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ABSTRACT 

We report a sensitive search for the HCN( J=2— >1) emission hnc towards SDSS J114816. 64+525150. 3 
(hereafter: J1148+5251) at 2=6.42 with the Very Large Array (VLA). HCN emission is a star formation 
indicator, tracing dense molecular hydrogen gas (n(H2) > lO^cm"^) within star-forming molecular 
clouds. No emission was detected in the deep interferometer maps of J1148-I-5251. We derive a 
limit for the HCN line luminosity of ^hcn ^ 10^ K kms~^pc^, corresponding to a HCN/CO 

luminosity ratio of L^(^j^/Lqq<0.13. This limit is consistent with a fraction of dense molecular gas 
in Jl 148-1-5251 within the range of nearby ultraluminous infrared galaxies (ULIRGs; median value: 



^hcn/-^CO- 



+0.05N 



and HCN-detected z>2 galaxies (O.lTlgQg). The relationship between L'^q^ 
and LpiR is considered to be a measure for the efficiency at which stars form out of dense gas. In the 
nearby universe, these quantities show a linear correlation, and thus, a practically constant average 
ratio. In J1148-I-5251, we find LpmJ L'j^^^>6600. This is significantly higher than the average ratios 
for normal nearby spiral galaxies {Lpm/ L'^fjj^=580^27n) ^^"^ ULIRGs (74015q^), but consistent with 
a rising trend as indicated by other z>2 galaxies (predominantly quasars; 1525tl75°)- It is unhkely 
that this rising trend can be accounted for by a contribution of active galactic nucleus (AGN) heating 
to LpiR alone, and may hint at a higher median gas density and/or elevated star-formation efficiency 
toward the more luminous high-redshift systems. There is marginal evidence that the Lfir/^hcn 
ratio in Jl 148-1-5251 may even exceed the rising trend set by other z>2 galaxies; however, only future 
facilities with very large collecting areas such as the Square Kilometre Array (SKA) will offer the 
sensitivity required to further investigate this question. 

Subject headings: galaxies: active, starburst, formation, high redshift — cosmology: observations — 
radio lines: galaxies 



1. INTRODUCTION 

High redshift galaxy populations arc now being de- 
tected back to 780 million years after the Big Bang (spec- 
troscopically confirmed: z=6.96; lye etal. 2006. ), prob- 
ing in to the epoch o f cosm ic reionization (e.g.. Fan etal. 
120061: Hu & Cowie|2006D. Many of these very distant 
galaxies show e videnc e for star formation activity (e.g., 
Taniguchi etal. |2005D. Some are even found to be hy- 
perlum inous infrar ed ga laxies (HLIRGs; Bertoldi et al. 
I2003at Wang et al. I2007D with far-infrared (FIR) lumi- 
nosities exceeding 10^'^ Lq, suggesting vigorous star for- 
mation and/or AGN activity. To probe the earliest stages 
of galaxy formation and the importance of AGN in this 
process, it is necessary to study the star formation char- 
acteristics of these galaxies. 

A good diagnostic to examine the star-forming envi- 
ronments in distant HLIRGs are observations of molecu- 
lar gas, the fuel for star formation. The by far brightest 
and most common indicator of molecular gas in galaxies 
is line emission from the rotational transitions of car- 
bon monoxide (CO), which was detected in ~40 galax- 
iesat high redshift (z>l; see Solomon & Vanden Bout 
|2005| for a review). These observations have revealed 
molecular gas reservoirs with masses of >10^'^Mq in 
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these galaxies, even in the highest re dshift quasa r known, 
J1148+5251 at z=6.42 (Walter et al. [200112001 Bertoldi 
etal. l2003bD . 

Although CO is a good tracer of the total amount of 
molecular gas in a galaxy, due to the relatively low crit- 
ical density of ~ 10^ — 10'^ cm^'^ required to coUi- 
sionally excite its lower J transitions, it is not a reliable 
tracer of the dense molecular cloud cores where the ac- 
tual star formation takes place. Recent studies of nearby 
actively star-forming galaxies have shown that hydro- 
gen cyanide (HCN) is a far better tracer of the dense 
(riHa ^ 10^ — 10^ cm""^) molecu lar gas where stars actu- 
ally form (e.g. Gao & Solomon l2004al l2004bl . hereafter: 
GS04a, GS04b). In the local universe it was found that 
the HCN luminosity (incN) scales linearly (unlike L'^q) 
with the FIR luminosity (I/fir) over 7-8 orders of magni- 
tude, rangi ng from Galactic dense cores to ULIRGs (Wu 
et al. [20051) . As -^^fir traces the massive star formation 
rate (unless AGN heating is significant), this implies that 
HCN is also a good tracer of star formation. 

HCN has now al so bee n detected in five galaxies at 
j>2 (S olomon et al. 120031: Vanden Bout etal 200J; Car- 



iUi et al. l2005l hereafter: C05; Wagg et al. [20051; Gao et al. 
[20?)7l hereafter: G07). Adding a number of upper limits 
obtained for other high-z galaxies, these observations in- 
dicate that the more luminous, higher redshift systems 
systematically deviate from the linear Ly^^j^-LpiR cor- 
relation found in the local universe (G07), and hint at 
a rising slope of the relation toward high Lfir and/or 
z. To further investigate this apparent non-linear, rising 
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Fig. 1.— VLA observations of HCN( J=2^1) (left) and 23.9 GHz 
continuum emission (right) towards J1148+5251 at a resolution of 
4.0"x3.3" (left) and 3.5"x3.3" (right; as indicated in the bot- 
tom left corner of both panels). The HCN{J=2— >1) line map is 
integrated over the full 25 MHz IF. No emission is detected in both 
maps. The cross indicates the geometrical center of the CO emis- 
sion. Contours are shown at (-3, -2, 2, 3)xcr (la = 15 /ijy beam~^ 
for the line map, and 14 ^Jy beam~^ for the continuum map). 

trend, our aim has been to extend the range of existing 
HCN observations beyond redshift 6 and to higher Xfir- 
In this letter, we report sensitive VLA"* observations 
of HCN(J=2^1) emission toward the z=6.42 quasar 
J1148-I-5251, the highest redshift source detected in CO. 
A previous, less sensitive search for HCN( J=2^1) emis- 
sion in this source has yielded no detection (COS). We use 
a concordance, flat ACDM cosmology throughout, with 
Ho=71 kms-iMpc-\ 17m=0.27, and f^A^O.TS (Spergel 
etal.[2007,). 

2. OBSERVATIONS 

We observed the HCN( J=2^1) transition line (z^rcst — 
177.2612230 GHz) toward J1148+5251 using the VLA in 
D configuration in 10 tracks on 2004 June 21 and 26 
(these observations were discussed by C05), and between 
2007 May 01 and 13. At the target redshift of 6.419, 
the line is shifted to 23.892873GHz (11.96mm). The 
total integration time amounts to 80 hr. Observations 
were performed in fast-switching mode using the nearby 
source 11534-^49311 (at 3.4° distance to J1148-h5251) 
for secondary amplitude and phase calibration. Obser- 
vations were carried out under very good weather condi- 
tions with 26 antennas [including 9 Expanded Very Large 
Array (EVLA) antennas in 2007] : the phase stability was 
excellent (typically <15° phase rms for the longest base- 
lines). The QSO 3C286 was observed as a primary flux 
calibrator. The resulting flux scale is accurate within the 
standard 15% uncertainty. 

In 2004, two 25 MHz wide intermediate frequency 
bands (IFs) with seven 3.125 MHz channels each 
were observed simultaneously, one centered at the 
HCN(J=2^1) line frequency (23.8929 GHz), and one 
centered at 24.0430GHz for continuum monitoring^. 
This leads to an effective bandwidth of 21.875 MHz for 
both the line and continuum observations, correspond- 
ing to 274 kms~^ at 23.9 GHz. For lines as broad as CO 

The Very Large Array is a facility of the National Radio As- 
tronomy Observatory, operated by Associated Universities, Inc., 
under cooperative agreement with the National Science Founda- 
tion. 

^ In fact, this is the center frequency of the HCO+{J=2— >1) line, 
which however was not detected. Neither was the continuum. 



TABLE 1 

Line luminosities in SDSS J1148+5251. 





[^Jy] 


L' 

[lO'^K kms-ipc2] 


Ref. 


HCN(J=2^1) 


(7 ± 15) 


<3.3 


1 


CO{J=1^0) 


<360 


<142 


2 


CO{J=3^2) 


570 ± 57 


26.4 ± 2.6 


3 


CO{J=6^5) 


2450 


26.9 ± 2.4 


2 


CO{J=7-»6) 


2140 


17.3 ± 2.4 


2 



Ref erenc es. — [1] This work, [2] Bcrtoldi ot al. I l2003bl) . [3] Walter 
otal. <2003f) . 



Note. — Bracketed number indicates nondetection. 

(279 kms-i FWHM; Bertoldi etal. l2003bl ). this setup 
would miss only ^15% flux in the line wings. In 2007, ob- 
servations were carried out in continuum mode, with one 
25 MHz wide IF centered at the HCN(J=2^1) line fre- 
quency (23.8929 GHz), and one 50 MHz wide IF centered 
at 23.7649 GHz to monitor the source's continuum^. 

For data reduction and analysis, the AXVS package 
was used. All data were mapped using the CLEAN algo- 
rithm and 'natural' weighting; this results in synthesized 
beams of 4.0 "x 3.3" for the line map and 3.5 "x 3.3" 
(~20kpc at z = 6.42) for the continuum map (see Fig.[T]). 
The final rms in the integrated HCN(J=2^1) line map 
is 15/iJy beam^^, and 14/iJy beam~^ in the continuum 
map. Averaging both uv datasets to a 'high-sensitivity' 
continuum map leads to an rms of 11 ^Jy beam~^. 

3. RESULTS 

No HCN( J=2^1) emission is detected in J1148+5251 
(Fig. m left). We derive a peak flux density of (7 ± 
15) ^Jy beam~^ at the source's CO position (see Ta- 
ble [J), setting an upper limit'' of 30 //Jy to the emission 
line peak flux. No continuum emission is detected at 
and/or close to the HCN(J=2^1) line frequency (Fig. [2 
right). From the 'high-sensitivity' continuum map, we 
derive a peak flux density of (11 ± ll)/xJy beam~^ at 
the source's position, setting an upper limit of 22/iJy 
to the 23.9 GHz continuum flux density. This is consis- 
tent with the model-predicted continuum level of <8 /iJy 
(Beelen et al.[2006). 

From our observations, we derive a limit to the 
HCN(J=2^1) fine luminosity of L'jjcn < 3.3 x 
10^ K kms^^pc^ (assuming a HCN/CO linewidth ratio 
of 0.67, i.e., the average value for the four z>2 HCN- 
detected quasars®; see Table [2] for references), corre- 
spond ing to 13% of the CO luminosity (Walter etal. 
I2003D . Using the FIR luminosity derived by Beelen et al. 
(|2006f) ^ we find Lfir/L'hcn>6600. 

4. ANALYSIS 

The nondetection of HCN(J=2^1) in J1148-f5251 
at the depth of our observations has several im- 

® This frequency avoids other potentially bright emission lines 
near the HCN(J=2^1) frequency, such as HCO+(J=2^1), 
HNC(J=2^1), and C2H(Af=2^1). 

^ In this letter, we quote 2a limits; however, note that all con- 
clusions hold for 3a limits. 

* This ratio lies between 57% and 80% for the z>2 sample, 
but is up to 100% in nearby starburst galaxies like NGC 253 (e.g., 
Knudsen etal. l2057h . 

^ FIR luminosities estimated from modeling the sparsely sam- 
pled FIR SEDs of high-2 galaxies are only accurate within a factor 
of 2. 
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TABLE 2 

Luminosity ratios at high z and average values. 



Source L-ei^jL'^^^ ■^'hcnZ-^'co 



SDSS J1148+5251 


>6600 


<0.13 


1,2,3 


Uloverleaf 


1305+360 


0.10+0.02 


4,5,6 


IRASF10214+4724 


2835+855 


0.18+0.05 


7,8,9 


VCV J1409+5628 


2615+860 


0.09+0.03 


9,10,3 


APM 08279+5255=^ 


1000 


0.27 


11,12,13 


SMM J16359+6612 


1550+550 


0.18+0.04 


14,15,16 


avg. 'normal' spiral (33) 
avg. LIRG (23) 


580_27o 


rid+° o^ 

'^•"^-0.02 
07+0 02 


17,14 
17,14 


avg. 2~0 ULIRG (9) 




n 17+0.05 
"■^'-0.08 


17,14 


avg. z~0 'air (33+23+9) 
avg. 2>2 galaxy (5) 


b80_3go 
1525lf,oo 


05+0 05 

n 17+0.09 

'^■^'-0.08 


17,14 
4-16 



galaxies; however, there is no indication that this fraction 
rises further toward the highest Lpm, or that it changes 
with redshift. 

4.2. 'Star Formation Law': The Lpm/L'^Q^ Ratio 

The FIR luminosity is thought to originate dominantly 
from dust-reprocessed light of young massive stars, and 
thus to be a good tracer of the star-formation rate (SFR) 
in a galaxy. Star formation in galaxies takes place in 
dense molecular clouds that are traced well by HCN 
emission, resulting in a linear correla tion b etween the 
FIR and HCN luminosities (Wu etal. I2005D . However, 
G07 have found a mild increase of L^m / L'^q^ with in- 
creasing FIR luminosity between normal spiral galaxies 
and ULIRGs [LFm/ L'^ci^=580tln (spirals) vs. 740t^°^ 
(ULIRGs) on average; see Table [2] . This trend appears 
to get stronger toward the even more FIR-luminous z>2 
systems {Lfiy{,/ L'^(^^=1525~!i\^^^ on average). The lower 
limit of Lfir/Lhcn>6600 obtained for J1148+5251 is 
consistent with such a rising trend toward high-z sys- 
tems, and extends it toward higher redshift (Fig. [21d). 
Due to the fact that the high-z sources also have system- 
atically higher ipiR than the nearby galaxies, it remains 
unclear whether we observe a rising trend with redshift 
or with LpiR (Fig. [St), or both. 

5. DISCUSSION 

5.1. Median Gas Density and Star Formation Efficiency 

Krumholz & Thompson ()2007[ ) argue that Lpi^/ L'^q^^ 
is expected to be higher for galaxies with a median 
(molecular) gas density n,„cd close to or higher than the 
critical density n^^f^ required for excitation of the ob- 
served HCN transition than for galaxies with lower rimed- 
In their case, they define star formation efficiency as the 
fraction of the mass that is converted into stars per dy- 
namical time of the system. Note that this is different 
than the star formation rate per unit total gas mass. 
They argue that the no n-linear relati on between LpiR 
and L/^jn (e.g., Kennicutt [l998al |1998^ GS04b; Riechers 
etal. [2006b) arises due to the fact that CO traces all gas. 
The star-formation rate is then dictated by the density 
n divided by the free-fall time (rg cx ri,~°'^), giving 
the standard Schmidt-law: star formation rate cx n^-^, 
or LpiR OC {L'qq)^-^. For molecules like HCN, which only 
trace the small fraction of dense gas clouds directly as- 
sociated with star formation in normal galaxies, Tff is 
roughly fixed by ridit- Hence the star formation rate 
shows a linear relationship with n, or LpiR cx (^hcn)^^- 
However, in extreme galaxies, where Timed m the molec- 
ular ISM approaches Tg again becomes relevant 
(i.e., HCN emission no longer selects just the rare, dense 
peaks whose density is fixed by n^™, but instead traces 
the bulk of the ISM, whose density can vary from galaxy 
to galaxy, and thus the variation of n and Tg re-enter 
the calculation), and the relationship approaches ipiR, 
cx (-^hcn)^'^ (and L'^qj^ cx L'^q)- Interestingly, current 
data show a marginal trend for a changing power-law in- 
dex at the highest luminosities of the type proposed by 
Krumholz & Thompson. This change in power-law index 
from 1 to 1.5 would suggest that, in these extreme lumi- 
nosity systems, nmcd approaches n^j^'^. More systems 
at high luminosity are required to confirm this trend of 
changing power-law index. 



Refere nces. — [1] This woric, [2] V^ alter etal. | |2003|), [3] Bee- 
lon etal. | |2006|). [4] S olomon etal. ||2003|1 . [5] Weifi etal. J2003r . [6] 
Riechers etal. l l2006aD . [7] Vanden Bout eta l J2004D . [8] D. Downes 
&: P. M. Solomon, in pr ep., [9 ] CariUi etal. I l2005f) . [101 Boelon etal. 
1 120041), [11] Wagg etal. Il20q 5|). [1 21 Riechers etal. (2006b), [13] Wcii3 



etal. I 2007 1. [14] Gao etal. II2007D , [15] K neib etal. (2005), [16] Knoib 



etal. I I2004I) . [17] Gao & Solomon I l2004bl) . 



Note. — Line luminosity unit is L^—K km pc . Bracketed num- 
bers are sample sizes ('all' corresponds to the full Gao & Solomon sam- 
ple). For the averaged samples, median values and Icr statistical errors 
of the cumulative distribution arc given. Where not quoted in the 
literature, 25% error arc assumed for LpiR. 

^ We adopt the extrapolated HCN(J— 1 — *0) luminosity from Weifi 
etal. (120071) . derived by assuming coUisional excitation only (uncer- 
tainties are model-dominated). Also, we only consider the 'cold, dense' 
gas and dust component in their model, which is assumed to give rise 
to the HCN emission. 

plications. For the following comparison with 
other galaxies (for which L^q and L'^qn 
given in the ground-state transitions), we assume 
that the HCN(J==2-^1) and CO(J==3-^2) lines are 
thermalized, so that iHCN(j=2-i)~^HCN(j=i-o) ^'^'^ 
-^co(J=3-2)~^co(J=i-o)- Considering the high CO ex- 
citation in this source (Bertoldi ct al. 2003b) and the fact 
that both lines arise from low J transitions, this is likely 
a valid assumption. 

4.1. 'Dense Gas Fraction': The -^^hcn/^co Roi^o 

CO(J=1^0) emission is considered a good tracer for 
the total amount of molecular gas in a galaxy, while 
HCN(J=1^0) emission is considered a good tracer for 
the dense peaks of the molecular mass distribution. The 
L'y^q^/L'qq ratio thus is considered a measure for the 
dense fraction of molecular gas in a galaxy^". On aver- 
age, this ratio is about four to five times higher in lo- 
cal ULIRGs compared to normal, nearby spiral galaxies 
[median value of O-Ht^^m (ULIRGs) vs. QMt^^*\ (spi- 
rals); see Table [2]. However, there does not appear to be 
a further increase toward the even more FIR-luminous 
systems at z>2 (0.17lQQg on average) within the sta- 
tistical uncertainties (all values rederived from GS04ab; 
GOT). The upper limit for L'^c^/L'^^ in J1148-h5251 is 
0.13, placing it below the median but within the range of 
values measured for ULIRGs and z>2 galaxies (Fig. [2^). 
Our observations thus confirm the finding of G07 that 
galaxies brighter than Lpir = 10^^Lq appear to have a 
higher fraction of dense molecular gas than normal spiral 

'^O We note that Gracia-Carpio et al. (|2006l ) have questioned the 
validity of ^'hcn/^CO ^ ^ tracer of the dense gas fraction. 
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Fig. 2. — Progression of the ^'hcnZ-^CO C'^'^iise gas fraction') and Lpi^/L^Qj^ ('star formation law') ratios witli redshift and LpiR. 
Data for nearby spirals (Lfir<1010-9 Lq i.e., Lir<10" Lq) and (U)LIRGs (Lfir>10^"-3 Lq; GS04ab; G07), and high-2 galaxies including 
J1148+5251 are shown (detections labeled, see Table |2]for references). In the left panel, down-arrows indicate upper limits for HCN, and 
up-arrows indicate that HCN was detected, but not mapped throughout the whole galaxy. In the middle and right panels, this is indicated 
vice versa. 

5.2. The Role of AGN Heating for Lfir 

Like most of the z>2 HCN-detected sources, 
Jl 148+5251 is a quasar. It has been found that, even for 
such strong AGN galaxies, the bulk of Lfir is likely dom- 
inatl y heat ed by star for mation in most case s (e.g. , Carilli 
et al. 120011: Omo nt et al. 120011: Beelen et al. 120061 : Riech- 
ers etal. I2006bf ). However, based on radiative transfer 
models of the dust SED of J1148-t-5251, Li etal. (2007!) 
argue that this source may currently undergo a 'quasar 
phase', in which AGN heating of the hot and warm dust 
contributes significantly to Lfir- If correct, this may be 
an alternative explanation for the elevated Lfir/Lhcn 
in this galaxy. The (rest-frame) IR properties (tracing 
emission from hot dust) of Jl 148+5251 are similar to 
those of other z>6 quasars with much lower Lfir (trac- 
ing emission from warm dust), and even to local quasars 
(Jiang etal. l2006f ). This supports the assumption that 
the hot dust in J1148+5251 is dominantly heated by the 
AGN; however, the lack of a correlation between Ljr 
and Lfir in quasars indicates that the warm dust may 
still be dominantly heated by star formation. Moreover, 
J1148+5251 follows the radi o-FIR correlation for star- 
forming galaxies (Carilli et al. |200"^ . which also suggests 
a starburst origin for the dominant fraction of Lfir- 

Furthermore, one of the z>2 HCN detections and some 
of the meaningful limits are submillimeter galaxies with- 
out a known luminous AGN, but are still offset from the 
local Lfir/LJjqj^ relation- It thus appears unlikely that 
AGN heating alone can account for the higher average 
Lfir/L^(jj.^ in the high-z galaxy sample. 

5-3. Implications for Future Studies 

Even when assuming the highest Lfir/L^^j^ of 2835 
found among all HCN-detected galaxies in Tabled the 
depth of our observations is sufficient to detect a galaxy 
with the redshift and Lfir of J1148+5251 (^=6.42) in 
HCN emission at a signal-to-noise ratio of >4.5. To first 
order, our lower limit thus is consistent with previous 
suggestions (G07) that Lfir/LJjqj^ ratios in high red- 
shift sources lie systematically above those for nearby 
galaxies. The scatter around this trend is still signifi- 



cant, and will primarily be improved by increasing the 
number of HCN-detected galaxies at high z. In addition, 
it will be important to improve on the main sources of 
error for the individual high-z detections (e.g., signal-to- 
noise limited HCN/ CO linewidth ratio, accuracy of the 
FIR SED fit, AGN bias of Lfir)- The statistical and in- 
dividual results, so far, would even be consistent with an 
even stronger increase in Lfir/L^qj,^ toward the high- 
est z and/or Lfir- Our study of J1148+5251 may hint 
at such an effect- Clearly, it is desirable to obtain more 
sensitive observations of this source to further investi- 
gate this issue- Due to its superior collecting area and 
high calibrational stability, the VLA is ideally suited for 
such a sensitive study- Although J1148+5251 is the most 
CO- and FIR-luminous z>& galaxy known, 80 hr of VLA 
observations were necessary to obtain the current limit - 
In a favourable case, the HCN(J=2^1) line may have a 
strength of about 1-5 times the current rms- To obtain 
a solid 5cr detection of such a line, of order 1000 hr of 
observations with the VLA would be required- Due to 
improved receivers and antenna performance, the fully 
operational EVLA will be by a factor of two more sen- 
sitive to spectral lines of several 100 kms""'^ width (such 
as in Jl 148+5251), but will still require long integration 
times- Studies of dense gas at z>6 thus appear to re- 
quire an order of magnitude increase in collecting area, 
such as offered by future facilities like the SKA phase 
I demonstrator (e.g., Carilli |2006[ ). which can serve as 
a low frequency counterpart to the Atacama Large Mil- 
limeter/submillimeter Array (ALMA)- 
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